VORTEX DYNAMICS FOR THE 
GINZBURG-LANDAU-SCHRODINGER EQUATION 

I> ' J. E. COLLIANDER AND R. L. JERRARD 

ON 

ov 

Abstract. The initial value problem for the Ginzburg-Landau-Schrodinger equa- 
tion is examined in the e — ► limit under two main assumptions on the initial 
data (f> € . The first assumption is that e exhibits m distinct vortices of degree ±1; 
i-^ ' these are described as points of concentration of the Jacobian [J<jf} of 4> e . Second, 

(vi . we assume energy bounds consistent with vortices at the points of concentration. 

Under these assumptions, we identify "vortex structures" in the e — > limit of (jf 
and show that these structures persist in the solution u e (t) of GLS e . We derive 
ordinary differential equations which govern the motion of the vortices in the e — > 
Q,, limit. The limiting system of ordinary differential equations is a Hamitonian flow 

■_! ' governed by the renormalized energy of Bcthuel, Brezis and Hclcin. Our argu- 

^3 . ments rely on results about the structural stability of vortices which are proved in 

a separate paper. 



> " 1. Introduction 

00 ' 

t*"" ■ We investigate the e — > behavior of the initial value problem 

CN ' 

rTq ftu e t -Au t + ^(\u e \ 2 -l)u t = 0, u e : T 2 x[0,T)i- 



GLSe \ u e (x,0) = ^(x), xeT 2 

exhibiting a finite number of vortices, under appropriate assumptions on the associ- 

^ ' ated Hamiltonian 

-<— > ■ 

* " (i.i) m = Jn U )dx; e\u) = i|Duf +^(i«f - 1) 2 - 

T 2 

The e — > asymptotics of minimizers of I e [u] subject to Dirichlet boundary con- 
straints were studied by Bethuel, Brezis and Helein in [[!]]. The corresponding para- 



bolic evolution problem was studied by Lin [|X3|| and Jerrard and Soner [llj| . These 
works identify "vortex structures" emerging in the e — ► limit whose location is gov- 
erned by a renormalized energy. The renormalized energy is obtained by removing 
the divergent part of P [u] . This note and the forthcoming paper Q establish similar 
results for the conservative evolution GLS f . 



The Landau theory of second order phase transitions [12| consists of expanding 
the energy in terms of a parameter which encodes the "order" in the phase and then 
exploiting energy properties to determine the evolution of the "order parameter". 
This theory was applied by Ginzburg and Pitaevskii and Pitaevskii [0 to ar- 
gue that the order parameter describing superfluid helium II evolves according to 
GLS € . In this context, P is the free energy and u € is the order parameter "which 
plays the role of 'the effective wave function' of the superfluid part of the liquid" 

l 
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||. The motion of u e under the GLS e evolution conserves 7 e [ii e ]. If we express 
u e (xi,x 2 ) = p(xi,X2)e t9( - Xl ' X2 \ with p, R-valued, then p 2 represents the density of 
the superfluid and DO is the velocity of the superfluid. Gross also derived GLS e 
as the Schrodinger equation for a wave function describing a system of interacting 
bosons. The equation GLS e is often called the Gross-Pitaevskii equation in the 
physics literature. 

Because we are interested in studying the motion of vortices in the e — > limit, 
we choose initial data (p e : T 2 i — ► C possessing vortices. This heuristic notion is 
implemented by imposing two conditions. First, we assume the measure [J4> e ] = 
det (D(j) e ) dx is close to a weighted sum of Dirac masses. More precisely, we assume 



;i-2) [J^] - nJ^dA 



i=i 



where d{ = ±1 with Yl^i = 0, and a, are distinct points in T 2 . 

8=1 

Second, we assume energy bounds. For example, certain of our results require that 

4> e satisfy 

(1.3) J e [0 £ ] <7rmlog- + 7 i. 

for some 71 G R. It follows from results that are cited and discussed in Section 3 that 
assuption ( |1.2|) implies the lower bound 

(1.4) P[(j) e } >7vm\og--C, 

where C is independent of e. Thus ( |1.3|) asserts that the energy of e is 0(1), given 
the prescribed configuration of vortices. Our strongest results need more stringent 
hypotheses, which amount to assuming that the energy is o(l), relative to the pre- 
scribed configuration of vortices. 

We briefly outline the contents of this paper. Since we have not yet defined many 
of the terms we are using, this description is necessarily a little impressionistic. 

Section 2 introduces notation and records some useful identities satisfied by solu- 
tions of GLS e . 

Section 3 states and briefly discusses some results providing conditions under which 
the Jacobian of a function u (which in our context may be interpreted as the vorticity) 
and the Ginzburg-Landau energy density concentrate around a collection of points. 
These results provide a key element in our analysis. The proofs are long and technical 
and will appear in ||. Section 3 also states some refined estimates related to the 
renormalized energy introduced by Bethuel, Brezis, and Helein fflj. These too are 
proven in Q. 

In Section 4 we prove that solutions of GLS e exhibit Lipschitz paths along which 
the energy and the vorticity concentrate. The main point is to control vortex mobility. 
We also describe weak limits of the "supercurrent" associated with solutions. 
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Finally, in Section 5 we give a complete description of limiting vortex dynamics, 
under somewhat stronger assumptions on the initial data. We also sharpen our 
characterization of the limiting behavior of solutions away from the vortices. 

We ultimately show that, under appropriate hypotheses, the limiting vortices solve 
a system of ODEs, which is exactly that satisfied by classical point vortices in an 
ideal two dimensional fluid. This result was predicted as early as 1965 by Fetter J/J, 
based on physical arguments, and more recently by Neu []IS]| and E ||, who arrived 
at the same conclusion through matched asymptotics. 

2. Notation and identities 

We introduce some vector notation and the quantity j(<p e )- We then define mea- 
sures [J0 e ] and p^ which will be used to locate the singularities of (j) e . Then a 
norm is introduced that permits us to say the singularities of two functions ifj € and 
(f) e are close if the associated measur es [Jip e ], p e ^ and [J4> e ], p e ^ are close. Next, 
we derive some identities for [Ju e ] and /i^ e assuming u e evolves according to GLS t . 
Finally, we introduce the renormalized energy W and a canonical harmonic map H 
with singularities of degrees {di, . . . , d m } located at the points {ai, . . . , a m } mapping 
T 2 \{a 1 ,...,a m }^S 1 . 

We introduce some notation for vector fields u : T 2 i — ► M. 2 mapping (x\, £2) G T 2 
to (u 1 (xi,X2), u 2 (x\,X2)) G M 2 . The scalar product in R 2 is denoted by "■", 

u ■ v = u v + u 2 v 2 . 

Sometimes we will view u G M 2 as u G C when we write iu or e ia u. These expressions 
are interpreted in the obvious way, e.g. 

' f 



iu = — Ju; J — , 



We define for u, v G 



p2 



j 



u x v = v}v 2 — u 2 v l = JijU % v\ 

V x u = d Xl u 2 — d X2 u x . 

Repeated indices are implicitly summed throughout this paper. Note that (iu) ■ v = 
u x v and (iu) ■ u = 0. For a scalar function (ft we define 

Vx = (cf) X2 ,-(j) Xl ). 

For a sufficiently differentiable u : T 2 1 — > IR 2 , we define 

j(u) = (u x u Xl ,u x u X2 ). 

We sometimes write j(u) = (iu) ■ Du. In the superfluid model, j(u e ) is interpreted 
as a current. If we write u = pe t6 , for p, 9 G K, then j(u) = p 2 D8. In particular, if 
\u\ = 1, j(u) is the phase gr adient. 
We define the signed Jacobian of u, 

(2.1) Ju = det Du. 
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The quantities j(u) and Ju are related by the identity 

(2.2) V x j(u) = 2u Xl x u X2 = 2Ju, 

so the signed Jacobian can also be interpreted as the vorticity. We write [Ju] to 
denote the distributional signed Jacobian, defined if j(u) G L l via 

(2.3) J<t>[M = lJ(Vx<l>)-j(u) : ^^(T 2 ). 

In particular, [Ju] is well defined for u G H l (M. 2 ). 
For a given u G H 1 ^ 2 ; R 2 ), define the measure 

VttA) = TT^—\ M«) dx; W{u) = \[Du[ 2 + -L(|u| 2 - l) 2 , 
| log e| J 2 Ae z 

A 

for subsets A <zT 2 . The renormalization factor i t , appears naturally upon consid- 
ering u = A smoothly cutoff in a ball of radius e centered at x = 0. 

We let M(T 2 ) denote the dual of C(T 2 ), i.e. the space of signed Radon measures 
on T 2 . Similarly, we let M X (T 2 ) denote the dual of C 1 (T 2 ), equipped with the dual 
norm. We also define a seminorm 

INIm1(T2) = swp < (j) dfj, : \\D(f)\\ Loo < 1, / = L 

If A*(T 2 ) = 0, we can compute H/^H^i^) by testing /x against functions such that 
j (/) dx = 0. In this case it follows that 

This is a consequence of the fact that ||</>|| c i( T 2) < C||.D</>|| Loo( - T 2n whenever f <fi = 0. 
In particular, if \x has the form 

n n 

(2-4) /i = ]T% - XX 

for some points £i, . . . , £ n , 771, . . . , r\ n G T 2 , not necessarily distinct, then Brezis, Coron 
and Lieb || show that 

II 

(2-5) II^IImW) = mi c n 5^l& - Vn(i)\ 

i=l 

where S n is the symmetric group of permutations of n objects. For measures of the 
form (|2.4j) we thus have 

n n 

c _1 rmnV|^ - r^ (i) | < ||//|| M i (T2) < min^|& - r^ (i) |. 
This estimate permits us to prove a fact illustrating the usefulness of the M 1 norm. 
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Proposition 1. Suppose that for every t G [0,T), /i t is a measure of the form 

n 

^2S ai M, for certain points ai(t), . . . , a n (t). The fj,r\ is a countiuous (resp. Lipschitz) 

1=1 

function from [0, T) into M 1 if and only if the points a.i{t) can be labelled in such a 
way th at ai(-) is continuous (resp. Lipschitz) for each i. 

Proof. For any s,t, the measure \i t — fi s has integral zero and so satisfies 



ir€S n 

i 



The lemma follows immediately. □ 

Remark 1. The M 1 seminorm can be interpreted as the minimum cost in a Monge- 
Kantorovich mass transfer problem. Indeed, when J dji = as above, « can be 
written in the form u = V\ — is 2 , where v\ and v 2 are positive, mu tually singular 
measures and ^i(T 2 ) = v 2 (T 2 ). Then ||/x||jgi is precisely the minimum cost of "trans- 
porting" V\ to z/ 2 , subject to an appropriate cost functional. See for example f or 
a more precise statement and more details. 

Remark 2. Any reasonable weak norm on measures would be equally suitable for our 
purposes. The M 1 norm is a convenient choice, but is certainly not the only possible 
choice. 

We employ O and o notation in some of the analysis below. We write, for example, 
O a ,6, c (l) to indicate a quantity is 0(1) with respect to the interesting limit e — > 0, 
with the implicit constant depending only upon the parameters a, b, c. 

We assume in this section, for fixed e > 0, that u solves 

iut — Am H — - \\u\ — l) u — 

and derive evolution equations for certain functions of u and its derivatives. These 
are mostly well known. We formally calculate the evolution equation for E e (u) = 

i|D.| 2 + ^(H 2 -i) 2 . 

d 1 

— E e (u) = u Xi ■ u Xtt + — (\u\ 2 - l)u-u t 



1 



= ( u Xi-Ut) x .-(u Xi:0i -—(\u\ -l))-U t 

so that 

(2-6) jW{u) = (u Xi ■ u t ) Xi 

using the equation and {iu t ) ■ u t = 0. 
Similarly, we calculate 

d l \ i 2 c \ A 

— -\u\ = u ■ Ut = —\iu) ■ Am 

dtT ' v ; 
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to obtain 

(2.7) dt2^ = ~^ iU ^ ' Ux ^ x > = _dlV J< ^' 
Next, we consider the evolution of a component of j(u). 

J t i k ( u ) = jj.(™) ■ u Xh = (iu t ) ■ u Xk + (iu) ■ u tXk 
= (iu t ) -u Xk -u- {iu tXk ) 

1 I ,2| |2 
= u x ] x ] ' Ux k — U ' UxjXjX k + ~2~\ u \ \ U \x k 

= [u x , ■ u x . — u ■ u x , x .) H — -ml ■ 

V X k J-J x k x ]J x - 2' 'X k 

We rewrite as 

rf 1 

(2.8) —j k (u) = 2(u x • «. .) — (Iw-b-I — u-u x . x .) H nlwL • 

Note that this shows J j2 j(u)dx is time- independent 
We calculate the evolution equation of Ju using (|2? 

-J^VX^-M^., 

which we write out explicitly, 

(2.9) j t Ju = $jk(u Xk -u Xl ) xiXj - 



Finally, we express (|2.9| ) in a weak form. Multiply Q2.9Q by r\ G C£°(T ) and 
integrate to find 

d 

It 2 Jt 2 

This is the key identity in our analysis of vortex dynamics of GLS e . 



(2.10) — J r] Ju dx = J r)xjxJjkU Xk ■ u Xl dx. 



Remark 3. The preceding calculations apply equally well to u e (-) G H l (T 2 ). For 
e > fixed, GLS e is a defocussing nonlinear Schrodinger equation. Bourgain has 
established [Q] global wellposedness for GLS e below H 1 . We validate the preceding 
calculations for u e G H 1 ^ 2 ) as follows using various aspects of Bourgain's result. 
By continuous dependence on the data, a different solution u e (t) is close to u e (t) in 
i^ x (T 2 ) provided the corresponding initial data <j) e and (fi e are close in if 1 (T 2 ). Let 
cf) e be a smooth approximator to e . The preceding calculations apply to u e since it 
remains smooth for all time. The various identities above, in particular ( |2.10| ), are 
then validated for u e G // 1 (T 2 ) by considering a sequence of smooth approximators. 

We reformulate for the torus T 2 some of the notions introduced by Bethuel, Brezis 
and Helein |l| for the study of Dirichlet minimizers of I e [u] . Let F solve 

(2.11) AF = 2tt(5 - 1) on T 2 . 
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with the normalization 

(2.12) lim(F(» -log|x|) = 

For points ai, ..., a m G T 2 (sometimes denoted by (a)) and nonzero integers d\, ..., d r , 
satisfying J^ o?j = 0, define the renormalized energy 

(2.13) W(a,d) :=-'K^d l djF{a i -a j ). 
Define 

m 

(2.14) <$> = ^diF{x - en) 



i=l 



and notice that 



m „ 



The following proposition is proved in [ED] (Theorem 1.1., Step 2 of proof). 
Proposition 2. There is a map H G C^JT 2 \ (a); S 1 ) n W^T 2 ; S 1 ) satisfying 

(2.15) j(i?) = -V x $. 

Moreover, if H is any other function satisfying ( |2.15| ), then H = e ia H for some 

a6R. 

It follows from ( |2.15| ) that H satisfies 

(2.16) divj(H) = 0, 

m 

(2.17) 2[JH] = V x j(H) = 2ixY,di5 av 

i=l 

in the sense of distributions and pointwise away from the singluarities (a). The map 
H also satisfies 

(2.18) / j{H)dx = 0. 

JT 2 

The map H has singularities of degree (of) located at the points (a) as seen in ( |2.17|) . 
The identity (|2.16| ) reveals that H is a harmonic map into S 1 since j(H) = D6 when 
we express H = e %e . Therefore, we refer to H as the canonical harmonic map^\ with 
singularities of degree (d) located at the points (a). We will sometimes highlight the 
dependence of H upon (d) and (a) by writing H (x; a, of). 

For a collection of points a\ , . . . , a m G T 2 define 

^^We have taken the name "canonical harmonic map" from the work of Bethuel, Brezis, and 
Helein on the Dirichlet problem, where it is entirely appropriate. In the periodic context it is clearly 
something of a misnomer, since H is not unique. Nonetheless, it seems easiest to use the familiar 
terminology. 
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We will be interested in the case 

(2.19) < p < -min |a» — a,j\. 

A configuration of vortices (a), (d) determines (up to a constant phase factor e ta ) a 
canonical harmonic map H(-, a, d). The map H has logarithmically divergent Dirich- 
let energy. The renormalized energy W(a, d) describes to leading orde r the finite part 
of the energy associated with the configuration (a), (d). The following proposition, 
essentially proven in (Theorem 1.7), makes these observations precise. 



Proposition 3. Let H = H(-,a,d). For p satisfying ( 12.191) , 



/ -\DH\ 2 dx = mnlog- + W{a,d) + 0(p). 
J T 2 p 2 p 

We define, following Bethuel, Brezis, and Helein |]J, the quantity 

(2.20) I(e,p) = min I / E e (u)dx : u G H x {B p ), u(x) = ^- for x e dB t 

[Jb p m 

It is shown in |l] that 

(2.21) J(e,p)=7rlog(^) + 0(l) 

as e — * 0, for p fixed. 

An easy adaptation of the construction given in ([[I]], Lemma VIII. 1) establishes 

Proposition 4. Given any distinct ai,...,o m G T 2 and di,...,d m G {±1}, there 
exists a family of functions v e G if 1 (T 2 ;M 2 ) such that 

[Jv e ] —*■ y^Tcdi5 ai weakly in M 

»=i 

and for every p > 

/ E e [v e ]dx < m{n log - + J(e, p)) + W{a, d) + Cp + o(l) 

Jt 2 P 

as e — »• 0. 

Remark 4. The condition [Jf € ] — *■ X) 71 "^^ weakly in M implies the energy lower 
bound 

f E e (u € )dx > m ( nlog- + I(e, p) ) + W(a,d) + Cp + o(l) 

JT2 V P / 

for some constant C and every p > satisfying fl2.19Q . This follows from Theorem 
2 below. Therefore, the functions v e of Proposition 4 are asymptotically energy 
minimizing, for the given configuration of vortices. 
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3. Concentration and structure 

This section identifies two hypotheses on a function u e which imply detailed struc- 
tural properties of u e as e — > 0. The first hypothesis is that [Ju e ] is somewhat close 
to a weighted sum of Dirac masses. The second is an energy upper bound consistent 
with the number of Dirac masses. Theorem 1 says that these hypotheses imply the 
measures [Ju e ] and p e uC concentrate to Dirac masses and that u e is uniformly regular 
away from the points of concentration. A more stringent energy upper bound on 
u e controls the current j(u e ) and D\u e \ as stated in Theorem 2. The proofs, which 



appear in |4j], rely heavily on techniques developed in UlQ| . Similar techniques have 
been used by Sandier fl7|. 



Theorem 1. (Global structure) Suppose thatu G if 1 (T 2 ;R 2 ), and that there exist 
points xi, ..., x m G T n , integers di, ..., d m G {±1}, and e < r := | min^ \xi — Xj\ such 
that 

m 

(3-1) || [Ju]-irJ2di6 Xi || M i(T2) < goo r ' 

and 

(3.2) / E e (u) dx < -Km log (-) + 71 

for some constant j-^. Then there exists points a« G B r / 2 (xi), % — 1, ...,m such that 

m / 1 \ 

(3-3) I K - ^zZ 8 ^ I U/i(T2) < 7l J- , 

in 

(3.4) \\[Ju] - 7r^d^ ai || M i (T 2) < o 7l (l), 

Moreover, for p fixed and < p < | ; 

(3.5) / £ e (u)cfo<0 P71 (l). 

(3.6) \\Du\\ L2(T2p) <O p ^(l), 

Finally for any p G [1, 2), 

(3.7) ||Dm|| Lp(T 2) < O mi (l) 
and 

(3-8) ||jHIIl,(t 2 )<O p , 7i (1) 

We have some remarks before commenting on the proof. 
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Remark 5. The estimate ( |3.1| ) says [Ju] is 0(l)-close to a weighted sum of Dirac 
masses located at the points X{. The energy upper bound (|3.2|) allows us to locate 
points a, near Xi and conclude [Ju] is o(l)-close to a weig hted sum of Dirac masses 
located at the points Oj. This observation is exploited in the proof of Theorem 3 
below. 

Remark 6. Theorem 1 implies several other useful estimates. For u satisfying the 
hypotheses of the theorem, we have by (|3.2|) that 

i 

(3.9) ||M 2 -l|| L2(T2) <Ce(log^ + l)". 

Interpolating with ( |3.7D gives for any p < oo 

(3.10) II m IIlp(t 2 ) — C v - 

For any fixed p with < p < |, the estimate ( |3.6| ) implies 

(3.11) II °'^'llz / 1 (T 2 ) — P' 



Remark 7. Related results for vector fields from MJ 1 i — > M. n generalizing the two 
dimensional results above appear in Q. 

We briefly describe the idea of the proof of Theorem 1 appearing in [[|. The 
theorem follows from a localized version with the hypotheses 

TV 

\\\Ju] — TrdS T . Lri/o i \\ < r 

/ E e (u)dx < 7rlog- +7i. 

JB r {xi) e 

The Jacobian condition implies u must exhibit nonzero degree in many subsets of B r . 
One way this condition can be satisfied, which is consistent with the assumed upper 
bound, is if u has a single isolated vortex of degree d near the center Xj. Careful lower 
bounds show that, if u deviates significantly from this description — for example, if u 
has more than one vortex, or if the energy or vorticity of u is spread out too much - 
then the total energy of u would be too large, contradicting the energy bound. The 
main technical point is thus lower bounds. These are proved by identifying some set 
which contains all concentrations of vorticity, and then covering it by balls in some 
optimal fashion. One can use this covering to relate the energy of the function to the 
distribution of vorticity. As noted above, this sort of argument has appeared before 
101 and in the work of Sandier Wj . 



in 



The next result provides additional information about u e under a stronger energy 
upper bound hypothesis. 

Theorem 2. Suppose that u e G H l is a sequence such that 

in 

(3.12) [J u<i ] ~^ n /_jdi5 ai weakly as measures, 

«=i 
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and that there exists some 72 > such that 

(3.13) / E e {u e )dx < m ( n log - + J(e, p) J + VP(a, d) + Cp + ^ 2 + o(l) 
as e — > ; /or every p > 0. T/ien 

(3.14) \imsnp\\^-j(u^ - j(H)\\ 2 L2{T2) < C l2 

and 

(3.15) limsupH D\u e \ 1 1 ^,2 ^2 \ < C72 

for every p > 0. Here H = H(-;a,d) is the canonical harmonic map and C is a 
universal constant. 

Remark 8 . Observe that the conclusion of the theorem is ridiculous if 72 is negative 
which implies the claim in Remark 4 . 

We conclude this section with a brief discussion of the m = case. Consider GLS e 
with initial data <p e = p e e ld with p e , 9gR. The solution is 

! (x,t)=^e i ^(^ a - 1 )']. 



u 



,2 



If \(p e ) — 1| > 0(e a ) for a < 2 then u € oscillates rapidly in the t variable and 
u e — ^ weakly in L 2 (dxdt) as e — > 0. The canonical harmonic map in this case is 
H(x) = e ldl for some constant 9\ G [0, 27r). The rapid temporal oscillations require 
\(p e ) — 1| = 0(e 2 ) to have u e — > H in any strong sense. 

The assumptions ( |3.2|) and m = imply \(p e ) — 1| = 0(e). Contrast this with 
assumption (|3.13| ) under the conditions 72 = and m = which yields |(p e ) — 1| = 



o(e). 

The quantity j(u) is insensitive to the temporal oscillations in u since j(u) = 
j(e ia u) for aGl. It is therefore natural to expect j(u) to retain more information 
in the e — > limit. This expectation is verified later in Theorem 3(iv). 

4. Vortex paths and convergence 

Theorem 1 describes a function u e (through the measures [Ju e ] and p e u ) provided 
that u e satisfies the hypotheses ( |3.1|) and ( |3.2|) . The theorem identifies Dirac masses 
as the e — > limits of [Ju e ], p e u . The hypotheses ( |3.1|) and (|3.2| ) are chosen to allow 
for a similar description of u e (t), the solution of GLS e with initial data ii e (0) = e , 
at each t in a time interval [0,T). In this section, we exploit the identity (|2.10|) and 
the results stated in Section 3 to show the Dirac masses in Theorem 1 move along 
Lipschitz paths provided u e solves GLS t . The paths are characterized as solutions of 
certain ordinary differential equations in the next section under stronger assumptions 
on the initial data. 

We assume the initial data e G ff 1 (T 2 ; M 2 ) satisfies the hypotheses (|1.2j ) and ( |1.3| ). 
So, there exist distinct points {«i, . . . , a m } C T 2 , integers {dx, . . . , d m } with d{ = ±1 
and Yl di = 0, fixed independently of e such that 
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in 



(4.1) [j^^^diSat 



i=\ 



(4.2) I E e ((j) e ) dx< vrmlog-+7i, 



e 

T 2 

We will also assume 
(4.3) / j^)dx = 0. 

In light of (|2.8|) , this implies J j(u £ (t))dx = for all time. This assumption is not 
essential but simplifies the exposition. There are modifications to the statements of 
Theorems 3 and 4 if J j(u e (t))dx ^ which appear in Q. 

For e > 0, Bourgain's recent work || shows GLS e is globally (in time) wellposed 
below H 1 . More precisely, there exists Sq < 1 such that for all s > Sq there is a 
uniquely defined continuous map 

H s i — >X C C([0,oo);H s ) 

taking 

where u e is the solution of GLS e . 

We can now state our first main theorem. 

Theorem 3. Let u e be the solution of GLS e with initial data (fi e satisfying ( |4.1| ) and 
( |4.2| ) and ( |4.3| ). TTien, a/ifer passing to a subsequence as e — > 0, there exists a T > 
(independent of e) and Lipschitz pat hs a^ : [0, T) \ — > T 2 , aj(0) = ctj such that 

m 

(i) [Ju e (t)] — *■ 7r^di<5 ai (i) weakly as measures for all t G [0,T). 
j=i 

fiij A^m ~~ *" ^X^M weakly as measures for all t G [0, T). 

i=l 

(in) \u € (t)\ 2 ->• 1 in L 2 (dx) /or all t G [0,T). 

(w,) j(w e ) = (i"u e ) • -Dw e — ^ j(H)) weakly in L p (dxdt) for all 1 < p < 2 where 
H(-,t) = H(-,a(t),d) is the canonical harmonic map with singularities of degree 
(d) located at (a(t)). 

Proof. 1. Define for r = min^ j\\ a i ~ a j\ with 0<e<r < 1 the quantity 

T £ = sup [t > : \\[Ju e (s)} - [J0 £ ]|| M1(T2) < ^r for all < s < t] . 
Recall that Ju = det Du = ul,ul„ — uiul. The estimates 



|[Ju e ]- [J01|| m i < ||Jw e - J(jf\\ L , < C\\u e - <j> e \\ Hl \\u e + 



\m 
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and continuity of the flow of (fi e i — ► u e {t) through H 1 guarantee T e > 0. 
2. Claim: For s,t satisfying < s,t < T e we have 

\\[J^(s)} - [Ju e (t)]\\ Mim < c\s -t\ + o 7l (l). 

Proof of Claim: The definition of T e guarantees for all t G [0, T e ) that the hypothesis 
(|3.1| ) of Theorem 1 holds if e is sufficiently small. Therefore, for each t G [0,T e ), we 



can find points a,(t) G Bz(ai), i — 1, . . . , m, for which 



(4.4) 



[Jli £ (t)] - TT^difi 



o,(t) 



8=1 



<o 7l (l), 



M!(T 2 ) 



by (|3.4j) . Of course, the aj(t) may depend upon e. So we can estimate 



< 



ir22di(5 a .( s) - 5 ai( t)) 



i=l 



\\[Ju e (s)] - [Jw\^j|| M i (T2) 

By ( |2~5|) , we can estimate by 

< 7r^|aj(s) - aj(t)| + o 7l (l). 



+ 0^(1). 



'11 



M^T 2 ) 



i=l 



The claim will be established once we show for i — 1, . . . , m, 

(4.5) |aj(s) - di(t)\ < c\s -t\ + o 7l (l). 

3. We prove ( [4.5|) by using the identity ( p.lO|) . Fix i and observe that aj(s), a^t) G 
Br(ai) for all s,t G [0,T e ). There exists an rj G C^°(S r (aj)) satisfying 

77(2;) = v ■ x for x G -Bar (fli(O)), vGS 1 

4 

and 

7r|a»(s) - Oi(t)| = vrdi / r/(5 ai(s) - ^(t))- 



The conditions on 77 guarantee that supp(D 2 r]) C -B r ( a i) \ B^on). Notice that 77 
depends upon the index z. 

Insert the function r\ described above into ( |2.10[ ) and again use ( [4.4| ) to observe 

t 

7r|oj(s) - Oi(*)| = Vxjxjjk u Xk ■ u Xl dx dr + o 7l (l). 

s B r (cti) 

The support properties of rj x . x , permit us to replace B r (a.i) by B r (cii) \ Bzr(a,i). 
Finally, we estimate by 

< \s-t\ \\D 2 r)\\ LOO sup \\Du(T)\\ L2(BAB3r) . 
re[s,t] t 

The size of ||Z) 2 ^|| LCX) depends upon r but is independent of e and ( |3.6| ) permits us to 
control the Du term by a constant independent of e, so ( |4.5| ) follows and the claim 
is proven. We also note that the claim implies T e may be taken independently of e, 
so we denote this quantity by T from now on. 
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4. The remaining convergence claims follow from the bounds stated in Theorem 1 
and passing to subsequences, except for (iii) which follows directly from (4.2). We 
prove (iv). Fix any p G [1,2). Since the conditions of Theorem 1 hold for every 
t G [0,T), we deduce from ( |3.8| ) that 

lliC u£ )llz,j>(T 2 x[o,T)) — Op,-n,TW- 
It follows, upon passing to a subsequence as e — > 0, that 

j(u e ) — *■ j weakly in L p (dxdt) 

for some j. We wish to identify j. 

Let (j) G C °°(T 2 x [0,T)). The identity (£77|) implies 

f M 2 

j(u e ) ■ Dtpdxdt = / <p t dxdt 

1 2 

0t _ (|w e | — l)dxdt 



as e — > for every t, by (iii). Therefore div j = 0. Moreover, from (i) we have 
V x j = 2 [J] <g> dt = 1-n J2 di8oi(t) ® dt weakly. 

Let H(x, t) = H(x, a(t), d). If we define V = j — j(H), we have 

div V = V x V = 

weakly. Let rf be a standard mollifier and set V s = V * r] s . The convolution here is 
in space and time. The above considerations imply 

div V s = V x V s = 

in T 2 for every t < T. Since V s is smooth, this implies V s (x,t) = g 5 {t). Letting 
5 — > 0, we find that V is also constant in x for each fixed t. For any fixed t we have 



V(x,t)dx = I (j -j(H))(x,t)dx 
j(x,t) dx 



lim / j(u € )(x,t)dx 



using flO|) . □ 

Remark 9. The proof given above may be iterated until 

T = inf{£ > : |Oj(t) — o,j(t)\ — > for some z 7^ j}. 
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5. Vortex dynamics 

Under the more restrictive upper bound on the energy ( p,13| ), we characterize the 
vortex paths Oj(-) as solutions of a system of ordinary differential equations that 
conserves the renormalized energy. More precisely, we have the folio wing result. 

Theorem 4. Suppose <j) e satisfies 

(5.1) [J^] ~^ ^y^di5 ai weakly as measures, 

i=l 

and for every p > as e — > 0, 



(5.2) 



/ E e (<j) € )dx < m U\og- + I(e, p) J + W(a, d) + Cp + o(l). 



Let {oi(-), . . . , a m (-)} and H be as in Theorem 3. Then the following statements hold: 
(i) For each i and for t e [0, Ti), 

'|ai = 2 J2 djV x F{a t - a 3 ) = -±diJD ai W(a, d) 

a*(0) = a*. 

i/ere Ti > can be taken to be the largest time such that the above ODE has a solution 

on[0,Ti). 

(ii) For every t e [0, Ti), 

r4i(0 - i(F) sirona/?/ m L? oc (T 2 \ (a(t))). 

(m_j For even/ p > and t e [0, Ti), 

lim min ||w e (-,£) — e ia /f (■, a(£), c()|| „ im , = 

Remark 10. For an example of a point vortex system which develops a singularity in 
finite time, see Marchioro and Pulvirenti [I3|. We expect, however, that for generic 
initial data the ODE does not develop singularities, in which case the above result is 
valid globally in time. 

Remark 11. The example given at the end of Section 3 shows that the sort of con- 
vergence in (Hi) is more or less optimal, unless (|5.2| ) is strengthened still further. 

Proof. 1. Let Oj(t), i = 1, . . . ,m denote the paths of Theorem 2. Let bi(t) denote 
the solution of the system 

(f t b i = 2J2d j VxF(b i -b j ) 
(5.3) <^ y.j& 

{ bi{0) = a t . 

A calculation shows that 

D ai W = -2i:^d l d j DF{a i - a,). 
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Therefore, the ODE in ( |5.3| ) may be reexpressed, as in (i), in Hamiltonian form 
showing that the renormalized energy W is conserved. This ODE system has a 
unique solution on a nontrivial time interval [0,T"). Let 

7\ = min(T, T) 

where T is as in Remark 9. Note that T\ is independent of e. We wish to show for 
all i that bi(t) coincides with ai(t) on the time interval [0,Ti). Observe that this will 
imply Ti = T' . 

For t e [0,71), let 

c(*)=X)i 6 *(*) -a *(*)i- 



It suffices to prove that, given any T < T\, we can find some small S(T) > and a 
constant C = C(T) such that 

d 



(5.4) 



dt 



at) < cm 



for a.e. t e [0,T] whenever £(£) < 5. We will show that ( p.4| ) holds at each point 
where Oj(-) is different iable for all i; by Rademacher's theorem, this condition is 
satisfied on a set of full measure. 

Fix T <T\. By Remark 9, there is some r = r(T) > for which 

(5.5) min \ai(t) — a>j(t)\ > 4r. 



A,t — Oi.t 



2. We use the fact that hi solves ( |5.3| ) and the triangle inequality to estimate 

x < Ei*. 



fit 



< 



2 E 



J^dyV x F(bi -bj)- ^djV x F(a % - o,-) 



+ E 



'j,t 



2^^-V x F(t»i- 
Term 1 + Term 2. 



We immediately dispose of Terra 1. Fix s < T and a pair of indices z 7^ j. Let 
/i = \(bi(s) — bj(s)) — (di(s) — a,j(s))\. Note that by assumption h < ((s) < 5. By 
Taylor's theorem, at the fixed time s, we have 



|Vx F(6 i -6 i )-Vx F( 



L> 2 F| 



r/,; — a_j) I < h max 

{a;:|a;— (a^— aj)|<h}' 

< C7C(s). 

The last inequality follows from ( |5.5| ) provided 5 < r. Therefore, Term 1 satisfies 

the desired estimate ( |5.4j ). 

3. We turn our attention to Term 2, which is a sum of terms (Term 2)j, with 



VORTEX DYNAMICS FOR THE GINZBURG-LANDAU-SCHRODINGER EQUATION 17 



i = l,...,m. Suppose that each function Oj(-) is different iable at s G [0,T). Fix 
rj G C^° such that suppijf) C B r ,fAai(s)) and 77(0;) = f • x in a neighborhood of eij(s). 
Here we take v G S* 1 to satisfy 



(5.6) 



(Term 2) ■ = d;z/ • a i<t 



(s)-2j2d J VxF(a l (s)-a J (s))). 



We will write e n for an appropriate subsequence as e — > 0. Since 

[Jw en (£)] — *■ 7rN^^i^ai(t) weakly as measures 
we can rewrite 

diV-diJs) = limdiV ■ -(ai(s + h) - a,i(s)) 
h-*o n 



lim lim — / (r]\Ju tn (s + h)] - T][Ju en (s)]) dx 

h—*0 n^oo7lhjj2 



s+h 



lim lim — - 

h— >0 n^oojch 



VxjxtJjkull ■ u^dxdt. 



We used (|2.10|) in the last step. 

Let H(x, t) = H(x, a(t), d). We reexpress the remaining term in 
1 which is stated and proven below, 



using Lemma 



d(u ■ 2Y^d,-V x F(a,i(s) — %(s)) 

s+h , >. 

lim- / diV- I ^djV x F(a,i(t) - afo)) J dt 



s+h 



l im J_ /" /" r, X:jXl J jk j k (H)j l (H)dxdt. 



Therefore, 



1 



s+h 



h— >0 n-tcojih 




(5.7) (Term 2), = lim lim — / / fc^ «» • t# - j k (H)j l (H)) dxdt. 



I 2 



Inside the integral, # = H(-, a(t), d) and Du e = Du € (t). 
On any set where \u\ > 0, we can decompose u x 



u 



■'■k 



iu 




iu 


u 




- — - 


■ u Xh 


- — - + 


- — - 


■ u Xh 


[\u\ 




\u\ 


\u\ 





u 

I u I 



(5i 



J (U) IU u 

+ M 



\u\ \u\ 



\x k 



u\ 
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Since \Du\ = a.e. on the set \u\ = 0, the representation above holds a.e. on supp(Du) 
and 

(5.9) K k -ui t = _LjV)iV) + \u%y\ xi . 

We will have proved ( |5.4j ) when we show 

s+h 

(5.10) lim lim - f f VxjXl S jk (|u e "L> e "U dxdt < C((s), 

s 

and 

s+h 

(5.11) lim lim \ [ [ VxjXl I jk (-^j k (u^)j l (u e ")-j k (H)j l (H)) dxdt < C((s). 

h^On-^oohj J T 2 \\u en \ J 

s 

These estimates will follow from the tight upper bound (|5.2|) (see also ( |3.13| )) on the 
energy and energy conservation. 

4. The renormalized energy W is conserved for solutions b(-) of ( |5.3| ). This is apparent 
from the fact that ( |5.3| ) can be written in Hamiltonian form. Also J T2 i? e (« € (-, t))(x)dx 

is conserved for solutions u e o i GLS e . Therefore, for every t < T and every p > 0, 
the upper bound (|5T2|) gives 



E e (u e (-,t))(x)dx = / E e ((f) t (-))(x)dx 

< m(n log - + I(e, p)) + W(a(0), d) + Cp + o(l) 
P 

= m(vr log - + I(e, p)) + W{b{t),d) + Cp + o(l). 
P 

Arguing as in the estimate of Term 1, we see that 

W(b(t)) - W(a(t)) < Cj2\Ut) ~ Oi(t)\ = C((t) 
provided 5 is small enough. Therefore 

(5.12) 

[ E e {u e {-, t)){x)dx < m ( vrlog - + /(e, p) ] + W{a{t), d) + Cp + C({t) + o(l), 

JT2 V P J 

as e — > for every p > 0. We have from Theorem 3 that 

m 

(5.13) [Ju En (t)] — 1 7r\^<ij5 ai ( t ) weakly as measures. 

i=\ 

The conditions ( |5.12| ), ( |5.13| ) are precisely the hypotheses of Theorem 3 with 72 

C((t). So, for every t G [s, s + /1], 

(5.14) limsupp|^"|(-,t)||^ (T?) < CC(t) 
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and 

(5.15) 



lim sup 



u L 



-j{u^)-j{H) 



< cat) 



L2(T2) 



These estimates allow us to prove (|5.10|) , (|5.11|) . 
We quickly estimate ( |5.10|) by observing 



(5.16) » < lim Jim \\D* V \\ L ^ 



s+h 



\D\u e 



\L*(T*) 



dt 



and using ( |5.14| ). 

5. We now establish ( p.ll| ). First, we show that 



(5.17) 



u' 



j(u e ") - j(H) -^ weakly in L 2 (Tj x [s, s + h}) 



To see this, note that r-^ij(u en ) is uniformly bounded in L 2 {T 2 x [s, s + h]) and hence 
converges weakly to some limit j. We know from Theorem 3(iv) that j(u e ") — *■ j(-ff) 
weakly in L p (dxdt) for all 1 < p < 2. We also know from Theorem 3 (hi) that 
\u e "\ rightarrow 1 strongly in L 2 (dxdt). Thus 

j(H) = weak L 1 lim j(u €n ) 



weak L 1 lim 



n— >oo \ \U 



3K U 



u 



( strong L 2 lim |-u en | ) I weak L 2 lim 

\ n—>oo J \ n— *oc 



jK 



+oo M e 



which proves ( 5.17| ). 

For fixed k, I, observe that the quadratic term in ( |5.11| ) can be written 



u c 



rj k (u^)-j k (H) 



M e 



-j\u^)+j\H) 



and reexpressed as 



f{u^)-j k {H) . 



Since ( j^nj k (u en ) — j k (H)) — ^ weakly in L 2 (T^ x [s,s + /i]) and j l {H) does not 

depend upon n, the second expression contributes nothing as n — > oo. The first 
expressi on is controlled using fl5.15| ). 

6. Since we have appropriately bounded Term 1 and Term 2, we have proven ( p.-l| ). 
Gronwall's inequality implies C, = which gives (i) of the Theorem. Since C = 0, 
( |5.15| ) implies (ii). We conclude by proving (iii). Fix t and p > 0. Let u e be a 



2(1 
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subsequence which converges in L 2 (T 2 ) to some limit u. We may assume, by ( ftp] 
that Du e — *■ D« weakly in L 2 (T 2 ). Also, (ii) gives 



liminf ||DM e ||i2( T 2) > lim 



e-^0 



e^O 



Jl" 



(/'- 



\Du\ 



L2(T2) 



L2(T2)- 



Therefore D« £ — > Du strongly in L 2 (T 2 ). Finally, since j(u) = j{H), Proposition 2 
implies u = e ia H for some aeK. 



□ 



We state and prove the lemmas employed above. 
Lemma 1. Suppose that rj e C 2 and that 

supp(r]) fl {oi, ..., a m } = {ai}; D 2 r/ = in a neighborhood of at 

Let H := H(-; a, d) be the canonical harmonic map. Then 



( Vxjx l S jk J k (H)j l (H) = diD-niai) ■ ( 2tt J" d 3 ; V x F 

h2 \ nii 



2 (J/t , 



Remark 12. This computation remains valid if di,...,d m assume arbitrary integer 
values, that is, if we lift the assumption that di = ±1 for all i. 

Proof. 1. We reexpress the integral in the lemma. Recall that j(H) = — V x $ where 
$ satisfies 



(5.18) 



A$ = ^27rd 4 <5 ai , 



»=i 



and, using (|2.14j) , we write 

(5.19) $(x) = d t F(x - a { ) + G(x) 
where 

G{x) = ^djFix-aj). 

Since j k (H) = Jt km ® Xm , we have 

(5.20) J jk j\H)j\H) = Sm$ Xn ® Xj . 
Fix any number p so small that D 2 7] = on B p (a,i). We have 

n X3Xl $ jk j k {H)j\H)dx = - / 7fc. ia .,Jj n $ SBn $ a 



T 2 



T 2 \B p (a,i) 



T 2 \B P 



»7x, hr&XnXjQxidx + / rix l h n ®x n Qx j v j dH 1 

dB p 
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where v = {y l ,v 2 ) is the outward unit normal to dB p . We recognize & XnX ®x — 
\{Q Xj <& Xj ) and integrate by parts again to find 

i 

.n jttI 



T 2 \B P 



/ ■n XlX J ln -$ Xj ® x .dx - I i] Xl J ln -$ < Xj $ Xj v n dH l 

dB p 

+ / T] Xl hn^x n ^x^ J dH l . 
dB p 

Since T] XlXn li n = 0, the integral over T 2 \ B p vanishes and we are left with two bound- 
ary integrals I p , II p . 
2. We calculate the boundary integrals. We begin with 

Ip = - / VxJln-^x^x^dH 1 . 
dBp 

By using ( |5.19| ), we observe 

(5.21) $ x . <$> Xj = F X] F x . + 2diF Xj G Xj + G x . G X] . 

Since F is even, the first term integrates to zero. We exploit the fact that G x . is nearly 
constant on B p (ai) to calculate the contribution to I p arising from the remaining two 
terms in ( |5.2i|) . The cross term contributes 

(5.22) 

-d l i lxi {a i )G Xj {a i )hn j F^dH 1 - *%, (o*) J,„ / ' F Xj [G Xj (x) - G Xj ( ai )]v n dH l . 

dBp dBp 

Since F ~ log \x — (h\, \F Xj \ ~ - on dB p and G is C 1 on B p , the second integral 
contributes O(p). The G x G x term contributes O(p) as well. 
Next, we calculate 

Up = / ^hn^x^x^dH 1 
dBp 

by expanding using ( |5.19| ). The F Xn F x . term again vanishes by symmetry. The 
G Xn G x . term contributes O(p) and the cross terms remain to be estimated. The first 
cross term gives 

(5.23) d^ia^G^a^hn j F^dH 1 + 0{p). 

dBp 

The second cross term contributes 

(5.24) di-q^GxSadhn J ' F^dH 1 + O(p). 



22 J. E. COLLIANDER AND R. L. JERRARD 



3. Since J gB F Xn vHE x = J gB F Xj v n dH l , the first term in |[22|) and (|Q3| ) cancel 



and the only remaining contribution is (|5.24|) . Finally, observe that 



(5.25) J F Xj v j dH 1 = fAFdx = 2n-2n 2 p 2 , 

dB p B p 



using (|2.11| ). Therefore 



(5.26) / Vx]Xl $ jk j k {H)] l (H)dx = 2nd i r ]xi {a i )S ln G Xn (a i ) + 0{p). 



T 2 

Since p can be taken arbitrarily small, we have proved the lemma. 

□ 
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